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We present numerical results obtained using a lattice-gas model with dynamical geometry defined by (Hasslachc: 
& Meyer 1998). The (irreversible) macroscopic behaviour of the geometry (size) of the lattice is discussed in 
terms of a simple scaling theory and obtained numerically. The emergence of irreversible behaviour from the 
reversible microscopic lattice-gas rules is discussed in terms of the constraint that the macroscopic evolution 
be reproducible. The average size of the lattice exhibits power law growth with exponent 1 /2 at late times. The 
deviation of the macroscopic behaviour from reproducibility for particular initial conditions ( "rogue states" ) 
is investigated as a function of system size. The number of such "rogue states" is observed to decrease with 
increasing system size. Two mean-field analyses of the macroscopic behaviour are also presented. 
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1. Introduction 

Lattice-gas automata (LGA) have successfully modelled a variety of systems including multiphase and 
multicomponent flow, chemical self-assembly and flow in porous media (Adler, et al. 1994, Rothman & 
Keller 1988, Love, et al. 2001, Boghosian, et al. 1996). Transformed into discrete quantum systems, suit- 
able generalisations of LGA have successfully modelled the Dirac and Schrodinger equations (Boghosian & 
Taylor 1998, Meyer 1996). All of these applications possess a fixed background geometry which is represented 
in the LGA models by a fixed lattice (and boundary conditions). In contrast, there are many natural systems 
with dynamical geometry. 

Only one lattice-gas model with dynamical geometry has so far been studied. Hasslacher and Meyer 
defined the only reversible lattice-gas automata with dynamical geometry in one dimension (Hasslacher & 
Meyer 1998). We are interested in modelling geometry intrinsically, rather than with reference to some am- 
bient space in which the system is embedded. Since it is the lattices in LGA which define the geometry, 
the model of Hasslacher and Meyer allows collisions in which lattice sites may be created and deleted. The 
restriction that the lattice-gas rules be reversible in timef restricts us to a single model in one dimension pro- 
vided we consider only two lattice vectors at each site. The classical, local, deterministic, single species LGA 
with two vectors per site and static geometry in one dimension is completely trivial as it can be interpreted 
as consisting of particles which simply move to the left or right without change even upon scattering. Never- 
theless, the one dimensional LGA with dynamical geometry is significantly more interesting (Hasslacher & 
Meyer 1998) . 

As it is the constraint of time reversibility which uniquely specifies the collision rules for the model we 
shall consider, it is worth considering carefully what time reversibility means in the context of the LGA. 
The evolution rule in all LGA has two phases. First each particle advects to the lattice point along its 
lattice vector. Second, the particles at each lattice point collide according to some local rule. The advection 
phase is trivially reversible and so it only remains to make the collision step reversible. In the deterministic 
models considered here this is achieved by making the collision a permutation amongst the states of the 
lattice site with the same values of the conserved quantities (in our case, mass and momentum). Physical 
time reversibility is not achieved by inverting all particle momenta and then running the forward evolution 
rule. Although advection and scattering are both reversible, these operations do not commute. Exact time 
reversibility is achieved by applying the inverse of propagation and collision in reverse order. 

f By which we will always mean local reversibility. 
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In the one dimensional dynamical geometry model the same advcct/collidc formalism for particle evolution 
is retained, but the collision rules allow for local changes of the lattice. An edge is created when a vertex 
splits into two vertices, and so the most local rule would depend only on the particle configuration at a single 
vertex. The exclusion principle restricts these configurations to consist of no particles, one particle, or two 
particles with opposite momenta. As explained in (Hasslacher & Meyer 1998) the constraint of reversibility 
allows only the possibility of splitting a vertex which is occupied by two particles, as shown in Figure 1. 
Reversibility then implies the dual edge deletion rule also shown in Figure 1. These rules are well defined 
as they cannot have overlapping domains of application. They define the simplest lattice-gas model with 
dynamical geometry. 

In keeping with the theme of this issue, we consider this model from the perspective of microscopic 
vs. macroscopic dynamics. When one seeks closed, effective macroscopic descriptions for systems out of 
equilibrium several problems arise. It is possible to define many macroscopic variables in terms of the plethora 
of microscopic degrees of freedom. It is often unclear which variables are necessary or sufficient to obtain 
a closed macroscopic description. We expect such macroscopic descriptions to be irreversible in general, by 
contradistinction with fundamental reversible microscopic laws. Finally, the domain of validity of macroscopic 
descriptions and when one may expect significant deviation from the "average" behaviour is often unclear. 
The model considered in this paper provides a simple example with which to consider some of these issues. 

The remainder of the present paper is organised as follows. In Section 2 we consider in very general 
terms the dynamics of the size of the lattice and introduce the idea of reproducibility of the non-equilibrium 
dynamics of the geometry to motivate the numerical studies of Section 3. We then analyse the model using 
two mean-field theory approaches in Section 4 before closing the paper with some conclusions and future 
directions. 



Figure 1. The scattering rules for our dynamical geometry LGA. Time runs up the page, and the lattice before and 
after the collision is depicted. Solid black dots indicate the lattice sites. The presence of a particle and its direction 
are indicated by an arrow. An edge is created when two particles advect to the same lattice point and an edge is 
deleted if after advection each endpoint is occupied by a single particle with outward momentum. Spectator particles 
are allowed in the empty vectors depicted here provided their inclusion does not violate the exclusion principle. 



Having defined our reversible lattice-gas model, we now wish to consider what its macroscopic behaviour will 
be. The "geometry" of the lattice in one dimension is trivial, the only macroscopic variable characterising 
the geometry is the size of the lattice. The time evolution of the size of the lattice will be dependent on the 
(reversible) microscopic evolution of the particle content. Denoting the size of the lattice and the microscopic 
state at time t by L(t) and s(t) respectively we may write the two evolution equations as: 



where £i(0) is the set of microstates corresponding to initial size Li and Ei(i) is the set of microstates 
resulting from t actions of the collision rule on each member of £^(0). Because the dynamics is deterministic 
the number of elements in £^(0) and Ej(i) is the same. We now define the property of exact reproducibility 
of the evolution of L(t). The set of microstates consistent with the new lattice size Lf(t) is E/(i), and exact 
reproducibility is the statement that all initial conditions with size Li(0) evolve to a final state with size 
Lf(t). In terms of the elements of the phase space this statement is: 




2. Macroscopic behaviour 



Si(0) 9 Si (0) ^ 8 f (t) e E<(t) 



(2.1) 
(2.2) 



Xi(t) C E f (t) 



(2.3) 
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For a lattice of size L containing p particles there are ( 2i ) initial conditions, and hence (2.3) implies: 




implying that if the evolution of the lattice size is exactly reproducible, then the evolution is irreversible: 

L(0) < L(t). (2.5) 

Clearly exact reproducibility is a very strong condition, and it is straightforward to construct microstates 
which violate this condition. We can, however, make a weaker statement based on the behaviour of a dis- 
tribution over the initial states. If we consider an initial distribution which is sharply peaked around -Lj(O), 
the weaker statement (which we will refer to as reproducibility) is that at later times the distribution will 
be sharply peaked around Lf(t). The argument above is repeated exactly, but now becomes a statement 
about the average lattice size over many independent realisations, and we must include the possibility that 
there will be "rogue states" in the tails of the distribution which lead to macroscopic behaviour far from 
the average. However, as the system size increases we expect the fraction of such microstates to decrease, 
until in the "thermodynamic limit" the probability of serious deviations from the macroscopic evolution is 
vanishingly small. This argument is presented and discussed in detail in (Garrett 1991). 

Having assumed that a simple macroscopic evolution rule for L(t) exists, and shown that its existence 
implies such a rule must be irreversible, we may now try and obtain the form of such a rule. In our numerical 
studies of the model we may therefore verify three things. Firstly, does the model geometry in fact possess 
a well defined reproducible behaviour? Secondly, if the mean behaviour is well defined, what form does it 
take? Thirdly, we expect there to be "rogue" microstates which will deviate from the macroscopic evolution. 
We can therefore explicitly count the number of such states as a function of increasing system size. 

3. Simulations 

Simulations were performed in order to obtain the average time evolution of the lattice size L(t). The 
system was initialised with 25 particles on a periodic lattice with 50 sites. The initial state of the system 
was chosen randomly, and 1000 different initial conditions were simulated for 10 6 time steps. The memory 
requirements for each realisation were sufficiently modest that individual simulations could be performed on 
a single processor. The large number of realisations required in order to obtain good statistics meant that 
the simulations were distributed across a 16-node Athlon Beowulf cluster. The lattice size was stored at each 
time step for each separate simulation. The mean lattice size as a function of time is shown in Figure 2. The 
mean lattice size as a function of time was fitted to a power law function given by Equation (3.1). 

L(t) = at b (3.1) 

The first 150000 timesteps were not included in the fit as this early time behaviour was not well described 
by a power law. The fits were obtained using the Levenberg-Marquardt algorithm and the errors on fit 
parameters were obtained by a careful ex post facto analysis of the behaviour of \ 2 near the minimum. The 
exponent obtained was 0.505 ± 0.02, for a fit with \ 2 — 0.47, consistent with a power law growth with 
exponent 1/2. 

The final numerical investigation we can perform is to search for the presence of "rogue" microstates and 
study their behaviour with increasing system size/ number of particles. In order to test the conjecture that 
the number of such rogue states will decrease with increasing lattice size, the time evolution of all microstates 
for initial lattice sizes 3, 4, 5, 6, 7, and 8 was computed for 10, 000 time steps. Two types of rogue states were 
found, states in which the lattice size remained constant, and states in which the lattice size oscillated. The 
fraction of rogue states (defined as those systems with lattice size less than 100 at time step 10000) was 
computed in each case and the results are plotted in Figure 3. These results are consistent with the number 
of rogue states becoming vanishingly small in the limit of large lattice size. 
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Time (time steps) 

Figure 2. Upper line: Ensemble average lattice size L(t). The lower line is a power law with exponent 0.5 and is 

included as a guide to the eye only. 
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Figure 3. Fraction of rogue microstates as a function of inverse initial lattice size for lattice sizes (from right to left) 

3,4,5,6,7,8. 
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4. Analysis 

(a) Mean-field theory - particles 
We consider one particular lattice site x. The probability that this site splits into two at the next timestep 

is: 

Pr(x splits at next timestep ) = Pr((x — 1, +) occupied and (x + 1, — ) occupied) (4-1) 

where (x,a) denotes a lattice site and momentum pair. If the density of right/left moving particles is p±, 
the "mean-field" approximation is: 

Pr(x splits at next timestep ) = p+p- (4-2) 

then the expected number of splits is E[+] = Lp + p^, where L is the number of lattice sites. Similarly, 
consider one particular edge e, between lattice sites x and x + 1. Then 



Pr{e deleted at next timestep ) = Pr((x, +) occupied and (x + 1, — ) occupied and (x — 1, +) 

not occupied and (x + 2, — ) not occupied) 



(4.3) 



Making the same mean-field approximation as above we obtain: 

Pr{e deleted at next timestep ) = p+p-(l — p+)(l — P-) (4-4) 
Combining the expected number of splits and contractions gives: 

E[AL] = E[+] E[-} 



= Lp+p-(p+ +p-- p+p-) 

writing p± = n±/L gives: 



(4.5) 



E[AL] = j^(Ln-p) (4.6) 
Where p = n + n- and n = n + + n_. Wc can therefore write a difference equation for E[L t \: 



E[L t+1 ] = E[L t ] + E[AL] 

P_ f ^ T i„ ^ ( 4 - 7 ) 
E[Lt 



= E[L t } + - E ^(E[L t ]n-p) 



We can convert this to a differential equation: 

dL p(nL — p) 
~dt = I? 

and integrate to obtain: 

(L-p/nf 3(L-p/n) 2 n i 3(L-p/n)n 2 



(4.8) 



2p p 



2 



n 3 

In (L - pin) = pnt + C (4.9) 

p 6 



This approach predicts an asymptotic growth of i 1//3 . This disagrees with our numerical data which shows 
an asymptotic growth of t 1 / 2 . The mean-field theory approach taken here is the natural thing to do if 
one is familiar with other types of lattice-gas model, where one defines average quantities which undergo 
macroscopic evolution. However, the lattice-gas model underlying our dynamical geometry model is non- 
interacting, and does not possess a local equilibrium. The average density for the model without dynamical 
geometry docs not posses a macroscopic evolution which can be obtained via the Chapman Enskog expansion. 
The validity of utilizing such mean quantities in the dynamical geometry case is therefore open to question. 

Secondly, the replacement of two-particle correlations with products of one-particle densities is appropri- 
ate if the collision frequency is rather high so that a particle rapidly becomes uncorrelated with neighbouring 
particles. In the dynamical geometry case, however, the density is a decreasing function of time and so this 
approximation, which underlies the mean-field theory, gets progressively worse at later times. 
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(b) Mean-field theory - lattice 



We have seen that a mean-field theory treatment based upon the average particle densities fails to 
reproduce the simulation results. As one goes further into the asymptotic regime the picture of an average 
density of particles becomes inappropriate because eventually the lattice is so large that the particles are 
very well separated. A treatment of the model in terms of the dynamics of the geometry, i. e. in terms of 
the particle separations, becomes appealing. 

Particle separations have two properties which determine their time evolution. They can be increasing, 
decreasing or constant, and they can be odd or even. From the point of view of the dynamics of the total 
lattice size, only decreasing separations are of interest. We see immediately that there is a natural timescale 
for the dynamics of separations. Immediately after a collision has occurred, the time to the next collision 
is given by half the size of the smallest decreasing separation. At the next collision the lattice will increase 
or decrease by one site. We introduce an ensemble with N members in which we label each realisation by 
i. If we denote the time to the next collision of the ith realisation by AT = Li/ct where Li is the size of 
realization i, and Cj > 1, then: 

Li(T + AT) = Li(T) + (N~ e - Nf°) (4.10) 

where Nf° (Nf e ) is the number of odd (even) decreasing separations of size Li in realisation i. Dividing by 
AT = Li/c, averaging over all members of the ensemble: 



E 



L(T + AT) - L(T) 
AT 



E 



(4.11) 



From our discussion of the connection between macroscopic reproducibility and irreversibility in Section 2 
we know that if the behaviour of the lattice size has a reproducible evolution, a necessary condition for the 
expectation value on the left hand side to have any meaning at all, the lattice size will be increasing in 
time on average. We therefore know that if the expectation value on the left hand side of equation 4.11 is a 
meaningful quantity then the right hand side will be positive. We now sum over i and introduce the number 
of increasing, decreasing and odd/even links in the ensemble via the numbers N AB , where A G {+, 0, — } and 
B G {o, e}. For N members of the ensemble and with P particles we have: 



NP 



E E NAB 

Ae{+,0,-} Be{o.e} 



(4.12) 



Using the above and replacing averages of products by products of averages we have 



E 



L(T + AT) -L(T) 
AT 



P 



E[a] (2N- e + N+° + N+ e + N 0e + N°°) E[a 



E\L] 



N 



EL] 



(4.13) 



As every quantity in the second term is positive, and the left hand side must be positive, the first term must 
be dominant here. We therefore neglect the second term and convert the above to a differential equation: 



(4.14) 



giving an asymptotic scaling of one half, consistent with the numerical data. 



5. Conclusions 

We have performed extensive simulations in order to characterise the late-time macroscopic behaviour of a 
one-dimensional lattice-gas model with dynamical geometry. We have adduced two mean-field arguments to 
explain the asymptotic growth. A mean-field theory based on the Boltzmann approximation for the particle 
densities yields an incorrect prediction of t 1 ' 3 growth. A similar argument applied to the separations of 
particles yields the correct exponent. 

As this model was originally motivated by the desire to describe interfaces in a two-dimensional multi- 
phase fluid intrinsically, the value of the growth exponent for the lattice size deserves some discussion. In 
binary immiscible fluids undergoing phase separation into droplet phases, three growth exponents have been 
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observed. An exponent of 1/3 is characteristic of a diffusive Lifshitz-Slyzov growth mechanism, and expo- 
nents of 1/2 and 2/3 have been observed in simulations including hydrodynamics. It is worth pointing out 
that two-dimensional lattice-gas simulations of droplet growth display a growth exponent of 1/2, where more 
mean-field oriented lattice-Boltzmann simulations obtained growth exponents of 1/3 (Wcig, ct al. 1997, Wag- 
ner & Yeomans 1998, Gonella, et al. 1999). This is at least an amusing coincidence, but we mention this 
mainly to caution the reader against reading too deeply into it. The physical interpretation of the lattice-gas 
particles in the one-dimensional model when considered as degrees of freedom of an embedded interface is 
unclear at best. 

This model possesses a generalisation to two dimensions which we are actively pursuing. As may be 
readily imagined, the two-dimensional model is considerably more complex. This is due in no small part to the 
presence of non-trivial local curvature. Work in progress includes characterising the hydrodynamic behaviour 
in a fixed, curved background geometry and producing efficient code for the simulation of dynamical geometry 
in two dimensions. 
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